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1 Introduction 



In the Standard Model it has been shown that bound states can be created 
for every quark but the top (see for instance fl], ||, ||, §, || and references 
therein). The latter possibility is ruled out due to the high value of the top 
quark mass, which is responsible for its short lifetime. The top quark decays 
directly into a W boson and a b quark before being able to create a bound 
state. In a recent paper || it has been shown that a finite probability of 
formation exists for a supersymmetric bound state made out of a squark 
stop and an antitop squark, for a certain range of t mass and some regions of 
the MSSM parameter space. In that, the possibility of signal detection at an 
e + e~ collider with LEP and NLC characteristics has also been investigated 
by means of a Breit-Wigner formula. 

A more refined result is needed in the threshold region which is charac- 
terised by low values of the squark velocity (3, i.e. 

/3 = yi--^«i . (i) 

For this purpose a by now standard Green function approach has been de- 
veloped (see |TJ, (J7|, H and references therein). 

We will assume that the supersymmetric bound state creation does not 
differ from the standard model CclS6 ; clS the relevant interaction is driven by 
QCD and controlled by the mass of the constituent squarks 0. For this 
reason the Schrodinger Green function technique is suitable for treating the 
problem of the scalar bound state. It will be used to compute the bound state 
effects on the cross section of stoponium near threshold. We will compare 
the results obtained in this manner to the Born cross section estimates for 
energies above threshold. The effects of the beam energy spread of the e + e~ 
colliders on the computed cross section will be discussed. 



2 The method 

The basic idea of the method is to consider the Schrodinger Green function 
equation |L| 

(H-£)G(x,y,£) = 5(x-y) , (2) 
where H is the Hamiltonian of the system 

H=-|| + nx) , (3) 
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and V(x) is the potential for the squarks. The imaginary part of the deriva- 
tive of the Green function given by @ taken at the origin is proportional 
to the cross section at threshold JTJ, [7|, H]. The finite width T of the state is 
taken into account by the substitution 



E -> E + iT 



Unlike the process 



the reaction 



e e 



tt 



e e 



tt 



(4) 

(5) 
(6) 



proceeds in P-wave, whereas top quarks are produced in S-wave configu- 
ration. This implies the Born level cross section to grow as a ~ (3 for top 
production (J5|), while for process (§) of scalars one obtains a slower rise, 
a ~ (3 3 . 

The threshold cross section of the process e + e~ — > tt is given by the 
following expression || pTOfl : 



a(e + e" -> tt) = R 



x 



q 2+ (^ + «e)Ql 



+ 



Mr 



s(s - w z 



sin2^ (s-M|) 2 + M|r| 



(7) 



where v e = —1 + 4 sin 2 6^, a e = —1, M 2 and Y z are the mass and the total 
width of the Z boson respectively. The charges are defined by Q 1 = —Q, 
Qz = (cos 2 9i~2Q sin 2 6y/) sin 20^ , with 6w being the standard electroweak 
mixing angle and 6~ t is the left-right mixing angle. 

The R term of equation (|7|) is obtained ]T|, 0, upon taking the imaginary 
part of the derivative at the origin of the Green function given by (0) 



R 



1536 tt 
81 m? 



x=0,y=0 

1 component of the 



(9* stands for the imaginary part, while G\ is the I 
Green function). 

For the purpose of our investigation, we will use a Coulombic potential 
for the Hamiltonian given in (Q) (here r = |x|) 



V 



4a, 
3 r 



(9) 



where a s is given by the QCD two-loop expression |TT 



2 



a s (Q 2 



Air 



(3 \og g 2 /A 2 



MS. 



1 



2/3 x log 


log 


\Q 2 / A hs}} 


Po log 


\Q 2 / A ws_ 





(10) 



and is calculated at a fixed value of the Bohr radius Tb given by 

3 



2mia s 



'IV, 



The validity of this choice has been shown (see || || HH an d references 
therein), and is essentially justified by the fact that the stop quark should 
be much heavier than all other quarks except (maybe) the top. The high 
value of mi implies from (|ll|) that the average distance between two squarks 
inside the scalar bound state is small, and therefore the dominant term of 
the potential is the Coulomb expression given in (H). This assumption for 
the bound state potential allows us to obtain an analytic expression for the 
Green function needed for the threshold cross section ||13|| . 

Following the above cited authors we introduce some standard notations: 
E = ^/s — 2m t ~ is the energy displacement from threshold, k 2 = —mfE, 
A = ?>a s mi/2 is the wavelength, and v = X/k is the wave number. Here the 
argument of a s is taken to be at the soft scale 1/rg . The finite width Y 
of the bound state is taken into account by means of the substitution given 
in (|). 

The expression for the I = 1 Green function for the Coulombic poten- 
tial (0) is given by 



367T 



A 2{k 2 - A 2 



2A 



+ In 



Gi(0, 0, k) 



2lE 



11 

IT 



+ 



(12) 



where 7^ is Euler constant (~ 0.57721), and ip\ is the digamma function, 
ipi(x) = d(\nT(x))/dx. The derivative at the origin of (|P^) is obtained by 
the simple multiplicative relation 



l^i Gl(x ' y ' fc) 



9Gi(0,0,fc) 



(13) 



x=0,y=0 



Some caveats (as described in []T3| ) have to be considered for the case 
of the P wave. There exists a constant linear term in the decay width F 
contributing to the I = 1 Green function that cannot be properly computed 
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in a purely nonrelativistic framework . The term independent from k that 
has to be added to the Green function ([12]) is given by 



0.185- ' 



36vr 



(14) 



and from |13[ the fif argument of the logarithmic term in (|T2"|) is given by 
0.13mf by means of an analysis of relativistic perturbation theory We note 
that the last two results have been obtained only for the top quark case, and 
we will assume that those results remain valid also for the stop quark which 
has a high mass, presumably close to that of the top quark. In any case 
the two aforementioned terms do not contribute much to the Green function 
estimates, since they are dominated by the leading k 3 term of (O). 



3 Results and discussion 

Our analysis of the threshold behaviour of the cross section should hold for 
a range of mass values and decay widths. For the mass range we will refer 



to the current stop mass value limits [[14]] and the LEP capabilities, while for 



the decay widths we have to take into account the formation requirements 
of the bound state ||. A criterion for the formation of bound states is that 
the creation of a hadron can occur only if the level splitting which depends 
upon the strength of the strong force between the (s) quarks and their relative 
distance 0, is larger than the natural width of the state. This means that, 
if 

A£ 2P _ip > r (is) 

where AE 2 p-ip = E 2 p — Ei P and Y is the width of the would-be bound state, 
then the bound state exists. Performing the analysis in this way allows us also 
to avoid dealing directly with several parameters of the MSSM model which 
are relevant to stop quark decay (apart from the stop mass), namely the 
ratio of the two Higgs vacuum expectation values tan /3, the Higgs-higgsino 



mass parameter /i and the wino mass M 2 fLOfl . The only constraint we have 
to impose is that the decay width be lower than about 1 GeV, necessary for 
the scalar bound state creation || regardless of the values assumed for the 
parameters mentioned above. 

As a first step, we check that the expression (0) used for our potential 
model @ is consistent with the Born cross section for E > 0. We see in 
figure ([I]) that in the non interacting limit with a s — > the cross section 
given by (|T2"D tends to the usual Born expression, and is zero for E < 0. This 



confirms the consistency of the expression found in (0). 
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In figures (QJ^) we present the results of the threshold cross section for 
a LEP mass range, mi = 60, 100 GeV. We have chosen the left-right 
mixing angle to be such that cos 2 ^ = 1. From (|7]) it is possible to see 
that the cross section minimum value is obtained for mi = 100 GeV at 
cos 6i ~ ±0.55, the point where the Z boson coupling vanishes, which does 
not much differ from the maximal value. As previously stated, assuming 
the decay width to be smaller than 1 GeV, we show for each chosen mass 
for widths of 10~ 3 , 10 -2 , 10 _1 and 1 GeV respectively. The cross sections 
are plotted against the threshold offset energy E, the centre of mass energy 
being given by the relation y/s = 2mi + E. 

For mi = 60 GeV, figure (fj) shows the structure of the discrete energy 
levels for E < versus decay width values given above. The maximal height 
of the peak, about 8400 fb, is obtained for the smallest width of 1 MeV. 
The shape of the peaks are similar to the one obtained by the Breit-Wigner 
formula, except for the resonance tails that are set higher with increasing 
energy. The height of the peaks decreases drastically as the binding energy 
reaches asymptotically the E = level. They tend also to accumulate and 
merge towards the E = value as the energy increases; this happens when 
the distance of the two resonance peaks is of the order of the decay width. 
For the Coulombic model the binding energy of the / = 1 level is given by 
the expression 

4 rniai , . 

9 n z 

and the resonance peaks merge when 



4 2 



1 1 

n 2 (n + 1) 



(17) 



This means that the the last visible peak has a quantum number n given by 

2n + 1 9r 



n 2 (n 2 + 1) 4m^a 2 



It 



we can see for instance in plot (|2|) that for T = 0.1 GeV the n = 3 peak is 
already barely noticeable. For T = 1 GeV , the limiting region of bound state 
formation ||, we see that there is no visible structure; even the first peak is 
smeared by the large width. 

The continuum energy region E > does not present any fine structure, 
and remains essentially the same for any given decay width. One important 
effect to be noted is the large difference of the Green function result with 
respect to the Born prediction for E > 0. This fact is shown in figure ([I]) 
where the two estimates of the cross section are compared. It is clearly 
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seen that the former is about one order of magnitude larger than the Born 
cross section. This can be understood by the fact that the Green function 
technique takes into account the interaction between the particles, and the 
contributions of the binding energies accumulate towards the E = energy 
level as described, thus affecting substantially the continuum region as well. 

In figure ([3D, we show the results for m t ~ = 100 GeV, with the same pa- 
rameters given in figure (0). The peaks are lower than the previous case, in 
particular the largest value of the cross section is obtained for T = 1 MeV 
at about 2400 fb. The position of the peaks are shifted towards lower val- 
ues of the centre of mass energy because the binding energy given by the 
expression ( |i~6|) is higher. 

For NLC energies, we present in figures and (^) results obtained for 
m f ~ = 200 GeV and mi = 500 GeV respectively. The parameters have been 
chosen to be the same as of the LEP case, and the qualitative behaviour of 
the threshold cross section is analogous. For mi = 200 GeV, the maximal 
value of the first peak obtained at T = 1 MeV is about 450 fb, while for 
mi = 500 GeV the highest peak reaches only 90 fb. The position of the peaks 
is shifted towards even lower energy values because of still higher values of 
bound state binding energy, as can be verified from ([[£]). 

We remark that all our results have been obtained for the Coulombic 
potential (|9|). It is known that in the threshold region there are singular 
Coulombic terms (a s //3) n which spoil the finite order perturbation theory. 
The resummation of these contributions have been done - see [O] and ref- 
erences therein - and they give small contributions and only modify slightly 
the Coulombic potential. The effect on the cross section is quite small, as 
can be seen from the plots in |Tj|, and does not change our estimates by 
more than a few percent. 

Another point concerns the validity of the Schrodinger Green function 
method (|T2|) . Since it is a nonrelativistic procedure, we have to ensure that 
the velocity of the squarks is low enough in order to make relativistic correc- 
tions negligible. From kinematical arguments, it is possible to give a bound 
for the maximal acceptable energy offset Emax ■ Assuming an upper value 
for the squark velocity, (3m ax , and the expression ([[]) together with the cen- 
tre of mass energy parametrisation, y/s = 2mj + E, one obtains by means of 
a series expansion in E 

Emax < mfiliAx ■ ( 19 ) 

In table ([!]) we present some estimates on Pmax and E M ax for different 
squark masses and some Lorentz 7 parameter values. We see that the limit 
of validity for the nonrelativistic equation (^) lies in a range of a few GeV 
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around the threshold, and naturally increases with larger squark masses; this 
can be intuitively understood since the heavier the squark is, the slower will 
it rotate inside the bound state 0, [|, ^. |(J . 

The results obtained so far have to be folded with the beam energy spread 
of the collider. For the LEP2 case, which has a beam energy spread of the 
order of 200 MeV |14|], even if the peak cross section is in the nb range for 
mi = 60 GeV @ the various resonance peaks are practically undetectable 
because their widths are much smaller than the typical beam energy spread 
(see also the discussion of The sole possibility of a width larger than the 
beam energy spread, T = 1 GeV already discussed, has no peaks as 

they have been smeared and thus no visible fine structure is envisaged. The 
situation does not change for the value of mi = 100 GeV; as we can see from 
figure (0) this situation is essentially the same as the previous one, and the 
peak cross section is even smaller than the former by a multiplicative factor 
of about 4. 

With the increase of the centre of mass energy (NLC case) the net result 
for the cross section detection is even worse than before. The beam energy 
spread is of the order of 6 GeV [15|] and, as seen clearly from figures (§) 
and ([5]), it is even larger than the energy range used for the plots by a factor 
of 4, thus making impossible the detection of any possible fine structure 
present in the threshold cross section. 



4 Conclusions 

In this letter, we have shown that the bound state effect on the threshold 
cross section of a scalar stop bound state is not negligible, at least for the 
case of a Coulombic potential. This effect turns out to be also dramatically 
different from the simple Born cross section results for E > [K], |j~6fl . 

However, because of the large beam energy spread of the present (LEP2) 
and future (NLC) e + e~ colliders, the possible structure of the cross section at 
threshold cannot be resolved. This confirmes our less refined Breit-Wigner 
approach of , and reinforces our previous result that the stoponium cannot 
be detected at the present and even future e + e~ colliders. 



Acknowledgements 

The author wishes to thank M. Antonelli for useful discussion and support 
during this work; and also O. Panella and P. Gensini for helpful hints. Many 
thanks also to Y. Srivastava for careful reading and patiently correcting the 



7 



manuscript. 



References 

[1] V.S. Fadin, V.A. Khoze, JETP Lett 46 525 (1987); 
V.S. Fadin, V.A. Khoze, Yad. Fiz. 48 487 (1988); 



[2] G. Pancheri, J.P. Revol and C. Rubbia, Phys. Lett. B 277 518 (1992). 

[3] J.H. Kiihn, E. Mirkes, Phys. Lett. B 296 425 (1992); 
J.H. Kiihn, E. Mirkes, Phys. Rev. D 48 179 (1993). 

[4] N. Fabiano, Eur. Phys. J. C 2 345 (1998). 

[5] N. Fabiano, A. Grau and G. Pancheri, Phys. Rev. D 50 3173 (1994); 
Nuovo Cimento A, Vol. 107, 2789, (1994). 

[6] M. Antonelli, N. Fabiano, Eur. Phys. J. C16 361 (2000). 

[7] V.S. Fadin, V.A. Khoze, Yad. Fiz. 93 1118 (1991); 

LI. Bigi, V.S. Fadin and V.A. Khoze, Nucl. Phys. B 377 461 (1992); 
LI. Bigi, F. Gabbiani and V.A. Khoze, Nucl. Phys. B 406 3 (1993); 
J.H. Kiihn, E. Mirkes and J. Steegborn, Z. Phys. C 57 615 (1993); 
W. Mddritsch, Nucl. Phys. B 475 507 (1996). 

[8] W. Kwong, Phys. Rev. D 43 1488 (1991); 

M.J. Strassler, M.E. Peskin, Phys. Rev. D 43 1500 (1991); 

M. Jezabek, J.H. Kiihn and T. Teubner, Z. Phys. C 56 653 (1992); 

Y. Sumino et al., Phys. Rev. D 47 56 (1993). 

[9] W. Beenakker, R. Hopker and P.M. Zerwas, Phys. Lett. B 349 463 
(1995). 

[10] K. Hikasa, M. Kobayashi, Phys. Rev. D 36 724 (1987). 

[11] W.A. Bardeen, A.J. Buras, D.W. Duke and T. Muta, Phys. Rev. D 18 
3998 (1978); 

W.J. Marciano, Phys. Rev. D 29 580 (1984). 

[12] N. Fabiano, Eur. Phys. J. C 2 345 (1998). 

[13] A. A. Penin, A.A. Pivovarov, Nucl. Phys. B 550 375 (1999); 
|hep-ph/ 9904278 . 



8 



[14] Review of Particle Properties, EuroPhys. Journ. C 3 1 (1998); 
[http: / / pdg.lbl.gov /| 

[15] Ed. R. Brinkmann, G. Materlik, J. Rossbach, A. Wagner, Conceptual 
Design of a 500 GeV e + e~ Linear Collider with Integrated X-Ray Laser 
Facility, Vol. 1 (1997); 
Ron Settles, private communication; 
Marcello Piccolo, private communication. 

[16] M. Drees, K. Hikasa, Phys. Lett. B 252 127 (1990). 



9 



Limits of nonrelativistic approach 

7 Pmax E MAX {60) E MAX {100) Emax{200) E max {500) 

1.01 0.140 1.18 1.97 3.94 9.85 

1.02 0.197 2.33 3.88 7.77 19.42 

1.03 0.240 3.44 5.74 11.48 28.70 



Table 1: Some estimates on (3m ax and Em ax as a function of the 7 param- 
eter for various stop masses, indicated in brackets, in GeV units. 




Figure 1: Comparison of the Green function method and the Born level 
expression for the cross section, in the limit a s — > 0. Here we assume that 
cos 2 61 = 1, mi — 100 GeV and a width of 0.5 GeV. 
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Figure 2: Cross section at threshold for various decay widths with mi = 
60 GeV, cos 2 61 = 1. The centre of mass energy is y/s = 120 GeV at 
threshold. 

Stoponium cross section at threshold 
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Figure 3: Cross section at threshold for various decay widths with mi = 
100 GeV, cos 2 9 1 = 1. The centre of mass energy is y/s = 200 GeV at 
threshold. 



11 



Stoponium cross section at threshold 

11^=200 GeV, cos 2 8, = 1 
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Figure 4: Cross section at threshold for various decay widths with mi = 
200 GeV, cos 2 Qi = 1. The centre of mass energy is ^Js = 400 GeV at 
threshold. 

Stoponium cross section at threshold 

m-=500 GeV, cos 2 a = 1 
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Figure 5: Cross section at threshold for various decay widths with mi = 
500 GeV, cos 2 9 1 = 1. The centre of mass energy is = 1000 GeV at 
threshold. 
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